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I.   I NTRODUCT I  ON 

The  understanding  of  the  kinematics  and  dynamics  of  the  impulsively 
started  flow  about  bluff  bodies  in  general  and  a  circular  cylinder  in 
particular  is  of  both  practical  and  fundamental  importance.  From  a 
theoretical  point  of  view,  analytical  and  experimental  results  enable 
one  to  understand  the  evolution  of  the  separation  of  flow  and  the  sym- 
metric vortex  development. 

From  a  practical  point  of  view,  the  understanding  of  the  character- 
istics of  impulsive  flow  enables  one  to  calculate  the  normal  and  side 
forces  acting  on  slender  bodies  of  revolution  at  large  angles  of  attack 
through  the  use  of  the  so-called  "impulsive  flow  analogy". 

The  present  investigation  addressed  itself  to  the  experimental 
determination  of  the  drag  (normal)  force  and  lift  (side)  force  acting 
on  smooth  cylinders  immersed  in  an  impulsively  started  flow  through  the 
use  of  a  vertical  water  tunnel. 


I  I.   PREVIOUS  STUDIES  AND  THE 
IMPULSIVE  FLOW  ANALOGY 

When  a  slender  body  is  inclined  to  a  flow,  it  experiences  a  compo- 
nent of  force  in  the  plane  in  which  the  inclination  occurs  and  a  com- 
ponent perpendicular  to  that  plane.  These  forces  will  be  referred  to 
hereafter  as  drag  and  lift  forces.  Earlier  works  were  mainly  concerned 
with  the  drag  forces  at  moderate  angles  of  attack.  The  flow  pattern  at 
small  to  moderate  angles  of  attack  consists  of  a  pair  of  axial ly 
oriented  symmetric  vortices  on  the  lee  side  of  the  lifting  body.  This 
physical  flow  field  was  approximated  by  the  so-called  NACA  vortex  model 
(See  Fig.  I).   Munk  PJ  was  one  of  the  first  to  study  the  problem  and 
to  introduce  the  cross-flow  approach,  which  has  been  the  basis  of  all 
drag-force  prediction  methods.   Munk's  work  was  primarily  concerned 
with  an  ideal  fluid  and  only  relevant  to  small  angles  of  inclination. 
Allen  and  Perkins  [2~]   considered  slender  bodies  at  large  angles  of  in- 
clination, where  the  viscous  effects  and  separation  are  important,  and 
modelled  the  cross-flow  development  by  analogy  with  the  flow  over  an 
impulsively  started  cylinder.   In  other  words,  they  have  hypothesized 
that  the  evolution  of  the  spiralling  vortex  sheets  with  distance  along 
a  slender  body  is  similar  to  the  evolution  of  symmetric  vortices  wi  th 
time  in  an  impulsively-started  flow  about  a  circular  cylinder.   This 
hypothesis  turned  out  to  be  a  very  fruitful  one  and  has  been  developed 
by  later  workers;  for  example  Sarpkaya  \JT\,    Kelly  [4]],  Hill  \J>~},    and 
Bryson  L~63  devised  various  approximate  discrete-vortex  models  to  study 
the  evolution  of  vortices  behind  an  impulsively  started  cylinder. 
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Figure  I.   Symmetric  vortex  separation  on  a  slender  body  of 
revolution  at  high  angle  of  attack,  NACA  model. 


11 


These  models  were  able  to  predict  in  an  approximate  way  only  the  early 
stages  of  motion.  These  investigators  were  not  concerned  with  the 
development  of  asymmetry  in  the  wake  and  subsequent  vortex  shedding. 

Vortices  peel  off  from  missiles  flying  at  high  angles  of  attack 
and  give  rise  to  significant  side  or  lift  forces.   In  spite  of  the 
realization  of  the  need  for  theoretical  analysis  and  experimental  data 
either  for  side  forces  acting  on  missiles  or  for  lift  forces  acting  on 
cylinders  in  impulsively  started  flow,  there  has  been  practically  no 
data  or  systematic  study. 

The  development  of  flow  about  a  cylinder  started  impulsively  from   ' 
rest  was  first  studied  by  Blasius  [7j.  He  has  shown  that  the  flow 
about  the  cylinder  does  not  begin  to  separate  until  a  distance  of 
s  =  0.351  c  is  covered  by  the  flow.   Blasius1  work  has  later  been  ex- 
tended by  Goldstein  and  Rosenhead  [8],  Gortler  [9]],  Schuh  DCG, 
Watson  (Jl!L  and  Wundt  C'2].  They  have  demonstrated  that  the  separa- 
tion depends  on  the  shape  of  the  body,  the  Reynolds  number,  and  the 
particular  time  dependence  of  flow  from  the  start  of  motion.   Follow- 
ing the  inception  of  separation  at  the  rear  stagnation  point,  the 
separation  points  move  rapidly  around  the  body  until  at  large  values 
of  time  they  reach  a  nearly  steady  state  position,  often  quite  close 
to  that  observed  experimentally  in  steady  flow  about  the  same  body. 

In  spite  of  its  relevance  to  the  impulsive-flow  analogy,  relative- 
ly little  experimental  work  has  been  carried  out  on  impulsively 
started  flow  about  cylinders.   Schwabe  [j3j  dragged  a  cylinder  of  9  cm 
diameter  in  a  shallow  channel  at  a  speed  of  0.8  cm/sec  through  the  use 
of  weights  and  pulleys  and  attained  Reynolds  numbers  of  about  600.   He 
developed  a  procedure  which  enabled  him  to  calculate  pressure 
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distribution  and  hence  the  drag  force  from  the  pictures  of  the  flow 
pattern.   Schwabe's  drag  coefficient  is  about  twice  as  large  as  the 
steady  state  value  and  is  still  increasing  after  the  cylinder  has 
displaced  a  distance  of  about  9  body  radii.  Approximate  numerical 
analysis  carried  out  by  Bryson  [_6]   shows  a  maximum  of  about  1.5  for 
the  drag  coefficient  at  the  time  when  the  cylinder  moved  about  4  body 
radii.  Subsequently,  the  drag  coefficient  drops  to  an  unreal i sti cal ly 
low  value  as  noted  by  Bryson. 

The  only  other  investigation  in  recent  times  was  undertaken  by 
Sarpkaya  Q3]  through  the  use  of  a  vertical  water  tunnel.  He  found 
that  the  drag  coefficient  reaches  a  value  of  about  1.5  at  a  relative 
fluid  displacement  of  about  8  body  radii.   In  spite  of  the  difficul- 
ties encountered  by  him  in  generating  an  impulsively-started  flow, 
his  results  have  clearly  shown  that  the  drag  coefficient  does  reach 
a  maximum  of  about  1.5  and  then  rapidly  decreases  to  a  steady  state 
value  of  about  1.2.  The  method  developed  by  him  has  been  further 
improved  and  utilized  in  the  present  investigation  not  only  for  the 
purpose  of  more  accurately  determining  the  exact  position  of  the 
maximum  drag  coefficient  but  also  the  evolution  of  the  lift  coeffi- 
cient. Consequently,  the  work  described  herein  presents  drag  force 
data  which  may  be  compared  with  those  obtained  previously  and  also 
lift  force  data  for  which  there  is  no  previous  work. 
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III.   EXPERIMENTAL  ARRANGEMENT 

A.  TEST  FACILITY 

The  experiments  described  herein  were  conducted  in  a  vertical 
water  tunnel  designed  by  Professor  T.  Sarpkaya  (See  Fig.  2).  A 
quick  release  valve  mechanism  which  opened  or  closed  with  the  help 
of  a  three-way  valve  was  mounted  at  the  bottom  of  the  tunnel.  When 
the  quick  release  valve  was  opened,  the  fluid  (water)  moved  down- 
ward through  the  tunnel  into  an  underground  reservoir.  The  partic- 
ular valve  design  is  shown  in  Fig.  3.   It  allows  an  initial  drop 
within  a  very  short  time  period  and  then  a  slow  further  opening. 
The  rate  of  initial  drop  as  well  as  that  of  the  subsequent  valve 
motion  is  controlled  by  a  piston  in  the  lower  part  of  the  valve  and 
the  viscosity  of  oil  in  the  liquid  chamber.   The  amount  of  the 
initial  drop,  the  area  of  the  opening  between  the  liquid  chamber  and 
the  upper  air  chamber  as  well  as  the  viscosity  of  the  oil  are  vari- 
able and  easily  adjusted  to  obtain  constant  velocities  at  desired 
rates.   Suffice  it  to  note  that  the  system  was  capable  of  producing 
an  impulsively  started  steady  flow  with  very  little  initial  disturb- 
ances at  desired  velocities.   In  fact  the  acceleration  of  the  flow 
was  confined  to  a  time  period  of  about  0.05  sec  for  steady  flow 
velocities  of  about  40  cm/sec. 

B.  ACCELERATION  AND  VELOCITY  MEASUREMENTS 

The  acceleration  of  the  fluid  was  monitored  by  means  of  a  differ- 
ential-pressure transducer  connected  to  two  pressure  taps  along  a 
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Figure  2.   Schematic  drawing  of  the  vertical  water  tunnel. 
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Figure  3.   Valve  assembly. 
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vertical  line  and  with  a  spacing  I.      The  acceleration  was  calculated 
from  the  fact  that  the  measured  differential  pressure  Ap  and  the 
acceleration  are  related  by 

P*  ^p  =  Ap  ( 1  ) 

in  which  p  is  the  density  of  water  and  —  represents  the  instantaneous 
acceleration.  Repeated  measurements  of  acceleration  have  shown  that 
the  velocity  reaches  a  constant  value  within  0.05  sec  and  the  accelera- 
tion remains  zero  thereafter. 

The  velocity  was  measured  through  the  use  of  a  variable  resistance 
probe.  The  change  in  resistance  of  an  about  150  cm  long  platinum 
wire,  placed  vertically  in  the  tunnel  away  from  the  tunnel  walls,  was 
recorded  by  an  amplifier  recorder  assembly  together  with  the  drag  or 
lift  force.  The  slope  of  the  elevation  versus  time  curve  yielded  the 
velocity  of  the  fluid  motion  through  the  use  of  the  calibration  data. 

C.   TEST  BODIES  AND  FORCE  MEASUREMENTS 

Two  smooth  aluminum  circular  cylinders,  60.6  cm  long  with  diameters 
62.23  mm  and  75.82  mm,  respectively,  were  used  in  the  experiments. 
The  cylinders  were  turned  on  a  lathe  from  aluminum  pipes.  The  length 
of  each  cylinder  was  such  that  it  allowed  a  gap  of  about  2  mm  between 
the  tunnel  wall  and  each  end  of  the  cylinder.  The  cylinder  was  pre- 
vented from  moving  towards  one  or  the  other  wall  by  means  of  small 
0-rings  attached  to  the  round  cantilever  end  of  the  force  transducers. 
A  double-ball  precision  bearing  was  inserted  at  each  end  of  the 
cylinder  in  aluminum  housings  which  sealed  the  cylinder  air  tight. 
The  outer  face  of  the  bearing  was  flush  with  the  cylinder. 
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Two  identical  force  transducers,  one  at  each  end  of  the  cylinder, 
were  used  to  measure  the  instantaneous  drag  and  lift  forces.  A  special 
housing  was  built  for  each  gage  so  that  the  gage  could  be  mounted  on 
the  tunnel  window  and  rotated  to  measure  either  the  drag  or  "the  lift 
force  alone.  The  bellows  protecting  the  strain  gages  were  filled  with 
Dow  Corning  3140-RTV  coating  for  water  proofing.  Then  the  ends  of  the 
bellows  were  sealed  air  tight. 

After  mounting  one  of  the  cylinders,  the  exact  angular  position  of 
each  gage  within  its  housing  was  determined  by  hanging  a  load  on  the 
cylinder  and  rotating  the  gage  until  the  lift  force  component  became 
zero.  The  final  positions  of  the  gages  were  marked  with  a  pin.   Finally, 
four  bolts  were  placed  on  the  gage  housing  to  hold  the  gage  rigidly 
in  position.  The  removal  of  these  bolts  and  the  set  pin  allowed  the 
rotation  of  the  gage  exactly  90  degrees  to  measure  the  drag  force. 
Thus,  both  gages  were  capable  of  measuring  either  one  or  the  other  com- 
ponent of  total  resistance.  Usually,  one  gage  was  set  to  measure  the 
drag  force  and  the  other,  the  lift  force. 

The  calibration  of  the  gages  was  accomplished  by  hanging  loads  at 
the  midlength  of  the  cylinder  after  setting  bofh  gages  to  sense  only 
the  drag  component.  The  calibration  was  also  checked  by  filling  the 
tunnel  with  water  and  recording  the  signal  generated  by  the  buoyant 
force.  This  procedure  provided  an  easy  check  on  calibration  and  was 
used  prior  to  each  experiment. 
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IV.   DISCUSSION  OF  RESULTS 

The  analog  data  for  both  the  lift  and  drag  forces  were  evaluated 
through  the  use  of  the  calibration  constants  and  plotted  in  terms  of 
the  normalized  distance  Ut/c  where  U  represents  the  velocity  of  the 
steady  flow;  t,  the  time  from  the  start  of  the  motion;  and  c,  the 
radius  of  the  cylinder.  The  drag  and  lift  coefficients  were  defined 
as  fol lows: 


r      .  Drag  Force  ,7. 

CD  "  1 —  1  U' 

j   p  L  D  U 

and 


Lift  Force  , ,. 

L  j   p  L  D  U 


in  which  p  represents  the  density  of  water;  L,  the  length  of  the 
cylinder;  and  D,  the  diameter  of  the  cylinder. 

A.   DRAG  COEFFICIENT 

Figures  4  and  8  show  the  drag  coefficient  as  a  function  of  Ut/c 
for  various  values  of  the  Reynolds  number  defined  by  Re  =  UD/v. 
These  figures  show,  regardless  of  the  Reynolds  number,  that  the  drag 
coefficient  always  rises  rapidly  to  a  value  of  about  1.55  at  Ut/c 
values  from  4  to  about  5  and  then  drops  rapidly.  The  rate  of  drop  as 
well  as  the  time  at  which  the  drag  coefficient  reaches  a  steady  state 
value  appears  to  depend  on  the  Reynolds  number.  At  relatively  lower 
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Reynolds  numbers  (Re  =  17,000),  the  drag  coefficient  drops  to  its 
steady  state  value  at  about  Ut/c  =  7  and  thereafter  fluctuates  about 
C_  =  I.I  as  a  result  of  the  lift-induced  drag  force  oscillations.  At 
relatively  higher  Reynolds  numbers  (Re  =  33,000),  the  maximum  peak 
in  the  drag  coefficient  occurs  at  Ut/c  =  5;  and  its  value  is  slightly 
lower  than  that  for  the  lower  Reynolds  numbers.  One  of  the  most  sig- 
nificant differences  between  the  evolution  of  the  drag  coefficient 
for  the  smaller  and  larger  Reynolds  numbers  is  in  the  manner  whereby 
the  drag  coefficient  reaches  its  steady  state  value.  As  noted  earlier, 
at  lower  Reynolds  numbers  Cn  rapidly  attains  its  steady  state  value 
at  about  Ut/c  =  7.  At  higher  Reynolds  numbers,  however,  Cn  first  de- 
creases from  its  peak  value  to  a  new  plateau  between  1.25  and  1.3  for 
Ut/c  values  from  7  to  1 3  and  then  decreases  once  more  to  its  steady 
state  value. 

The  foregoing  discussion  covered  only  the  general  characteristics 
of  the  evolution  of  the  drag  coefficient  with  time  without  making  an 
attempt  to  explain  the  reasons  giving  rise  to  them.  The  present  ex- 
periments as  wel I  as  those  previously  conducted  by  Sarpkaya  C33  show 
that  during  the  initial  and  symmetric  development  of  vortices  (as 
evidenced  by  the  absence  of  lift  force  as  will  be  discussed  later) 
the  vorticity  rapidly  accumulates  in  the  wake.  The  symmetric  vortices 
grow  to  sizes  considerably  larger  than  those  which  would  be  found  in 
the  later  stages  of  the  motion  where  the  vortices  shed  a  I ternat ing ly. 
The  accumulation  of  vorticity  reduces  the  back  pressure  and  gives  rise 
to  a  drag  coefficient  of  about  40*  larger  than  its  steady  state  value. 
In  other  words,  the  initial  rise  in  the  drag  coefficient  is  a  conse- 
quence of  the  symmetry  of  the  vortices  and  the  accumulation  of 


of  vortlcity  during  this  symmetric  growth.  Once  the  vortices  become 
asymmetrical  under  the  influence  of  ever  present  small  disturbances, 
one  vortex  moves  away  from  the  cylinder  causjng  a  rapid  drop  in  the 
drag  force  and  a  gradual  increase  in  the  lift  force.  Evidently  the 
rate  of  rise  and  the  peak  value  of  the  drag  coefficient  are  not 
materially  affected  by  the  disturbances  since  the  vortices  continue 
to  retain  their  symmetry  for  Ut/c  values  smaller  than  about  5.   It  is 
possible  to  show  that  the  vortices  must  acquire  certain  strength  and 
position  relative  to  the  cylinder  in  order  to  reach  a  state  at  which 
they  are  most  susceptible  to  disturbances.  Once  the  vortices  begin  to 
shed,  the  drag  coefficient  reaches  a  value  of  about  I.I,  a  value  which 
is  commonly  accepted  at  the  Reynolds  numbers  encountered  herein. 

The  peak  value  of  the  drag  coefficient,  its  time  of  occurrence, 
and  subsequent  decrease  to  steady  state  values  are  also  related  to  the 
character  of  the  fluid  motion  (laminar  or  turbulent),  and  hence  on  the 
dissipation  of  circulation.   It  is  a  well  known  fact  that  the  fluid 
motion  in  the  near  wake  is  turbulent  for  Reynolds  numbers  larger  than 
about  200*  The  rate  of  decrease  of  circulation  in  each  vortex  depends 
on  the  intensity  of  turbulence  and  hence  on  the  Reynolds  number. 
Numerical  analysis  through  the  use  of  discrete  vortices  C3j  has  shown 
that  the  larger  the  dissipation,  the  smaller  is  the  peak  value  of  Cn 
and  the  more  gradual  is  the  decrease  of  Cn  to  its  steady  state  value. 
Evidently,  a  vortex  pair  of  smaller  circulation  gives  rise  to  a  smaller 
drag  coefficient  since  Cn  is  proportional  to  the  square  of  T/Uc.   The 
experimental  results  reported  herein  are  in  conformity  with  the  fore- 
going conjectures.  At  higher  Reynolds  numbers  (Re  =  33,000)  the  peak 
value  of  C   is  about  1.5  and  occurs  at  a  Ut/c  value  closer  to  5,  in 
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constrast  to  a  CD-value  of  1.55  at  Ut/c  s  4  for  Reynolds  numbers  of 
about  16,000.  The  dissipation  of  vorticity  evidently  requires  that  the 
symmetric  vortices  be  fed  a  relatively  longer  time  period  CA(Ut/c)  =:  G 
in  order  to  allow  the  vortices  to  reach  their  maximum  strength.  As 
noted  earlier  the  vortices  gradually  become  asymmetrical  and  a  turbulent 
vortex,  being  inherently  more  stable  to  disturbances  relative  to  the  one 
at  lower  Reynolds  numbers,  lingers  somewhat  longer  in  the  vicinity  of 
the  cylinder.  These  ideas  have  previously  been  substantiated  through 
measurements  and  flow  visualization  by  Sarpkaya  L~3D.  The  extended  stay 
of  the  vortices  in  the  vicinity  of  the  cylinder  does  not  allow  the  drag 
force  to  drop  rapidly  to  its  steady  state  value.  Suffice  it  to  note 
that  the  evolution  of  the  drag  force  with  time  depends  on  the  Reynolds 
number.  The  larger  the  Reynolds  number,  the  larger  is  the  dissipation 
of  vorticity,  and  the  smaller  the  peak  value  of  CR.   It  has  been  noted 
previously  that  there  is  a  plateau  in  CR  for  Ut/c  values  between  7  and 
13  for  larger  Reynolds  numbers.   It  is  now  clear  that  that  plateau  is 
a  consequence  of  the  dissipation  of  circulation  in  the  vortices.   It  is 
also  clear  that  one  must  further  substantiate  these  ideas  by  carrying 
out  experiments  at  even  higher  Reynolds  numbers. 

B.   LIFT  COEFFICIENT 

The  lift  coefficient  versus  Ut/c  plots  are  presented  in  Figs.  9 
through  12  for  various  values  of  the  Reynolds  number.  The  Reynolds 
numbers  in  these  plots  do  not  exactly  correspond  to  those  in  the  drag 
coefficient  plots  because  of  the  fact  that  the  lift  and  drag  force 
runs  were  carried  out  independently.   It  was  not  possible  to  record 
simultaneously  the  elevation  and  the  lift  and  drag  forces  because  of 
the  lack  of  a  high-speed  three-channel  recorder. 
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A  quick  examination  of  the  lift  plots  show  that  the  lift  force  re- 
mains practically  zero  until  the  vortices  become  asymmetrical  at  Ut/c 
values  of  about  6,  the  amplitude  of  the  lift  force  varies  significantly 
from  one  run  to  another  even  for  identical  Reynolds  numbers,  and  the 
first  few  peaks  of  the  alternating  lift  force  do  not  occur  at  the  same 
Ut/c  values.   In  general,  however,  the  lift  coefficient  reaches  a  maxi- 
mum value  of  about  0.4.  This  value  is  in  conformity  with  those  obtained 
previously  in  steady  flows  at  the  corresponding  Reynolds  numbers  C 1 43- 

There  are  various  reasons  for  the  observed  behavior  of  the  lift 
force.  Firstly,  even  in  steady  flows  the  lift  coefficient  shows  large 
variations  at  a  given  Reynolds  number  for  all  Reynolds  numbers.  This 
is  primarily  because  of  the  inherent  instability  of  the  vortices,  the 
surface  condition  and  length-to-diameter  ratio  of  the  cylinder,  and 
the  lack  of  spanwise  coherence  of  the  vortex  cores.   Secondly,  in  im- 
pulsive flow  the  development  of  lift  is  strongly  dependent  on  the 
initial  disturbances  which  give  rise  to  asymmetry  and  subsequent  vortex 
shedding.  Consequently,  it  is  not  too  surprising  to  find  that  in  some 
runs  the  vortices  continue  to  grow  symmetrically  and  thus  result  in 
very  small  lift  forces.   In  some  runs  asymmetry  develops  rapidly  and 
the  lift  force  not  only  starts  at  smaller  values  of  Ut/c  but  also 
reaches  larger  amplitudes.   Since  there  is  no  way  to  control  the  magni- 
tude of  the  disturbances,  other  than  assuring  that  the  fluid  in  the 
tunnel  is  as  quiescent  as  possible  prior  to  each  run,  one  must  antici- 
pate relatively  large  variations  both  in  the  magnitude  and  the  time 
of  occurrence  of  the  peak  values  of  lift.   Suffice  it  to  note  that  the 
lift  force  is  a  fairly  random  quantity  and  that  randomness  is  certainly 
reflected  in  the  data  presented  herein.  The  data  further  helps  to 
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explain  the  random  nature  of  the  side  forces  measured  on  missile 
models  mounted  in  wind  tunnels.   It  appears  that  attempts  to  obtain 
a  unique  lift-force  variation  will  not  be  fruitful  in  view  of  the 
stated  randomness  of  the  lift  forces. 
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V.   CONCLUSIONS 

The  experimental  investigation  of  the  impulsively-started  flow 
about  smooth  circular  cylinders  in  a  vertical  water  tunnel  warranted 
the  following  conclusions: 

1.  The  drag  coefficient  reaches  a  peak  value  of  1.5  to  1.55, 
depending  on  the  Reynolds  number,  at  Ut/c  values  from  about  4  to  5; 

2.  the  drag  coefficient  eventually  drops  to  its  steady  state 
value  for  Ut/c  larger  than  about  7; 

3.  the  lift  coefficient  begins  to  develop  at  Ut/c  values  larger 
than  about  6  and  alternates  about  a  zero  mean  value  with  an  amplitude 
of  about  0.4.  The  magnitude  and  the  time  of  occurrence  of  the  peak 
values  are  strongly  dependent  on  the  magnitude  and  distribution  of 
the  initial  disturbances; 

4.  the  lift  and  drag  coefficients  both  appear  to  depend  on  the 
Reynolds  number.  Only  additional  experiments  at  higher  Reynolds 
numbers  could  determine  the  extent  and  form  of  the  dependence  on  the 
Reynolds  number. 
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